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ABSTRACT 


The past twenty years have seen an explosion in the realm of infrared deteetion 
teehnology fueled by improvements in III-V semieonduetor teehnology and by new 
semieonduetor growth methods. One of the fastest growing areas of this researeh 
involves the use of bandgap engineering in order to ereate artifieial quantum wells for use 
in Quantum Well Infrared Photodeteetors (QWIPs). QWIPs have an advantage over 
other infrared deteetors such as Mercury Cadmium Telluride (MCT) because they have 
larger bandgaps and are therefore stronger and cheaper to manufacture. This thesis 
introduces one method of “multi-color” detection through the use of an asymmetric 
quantum well structure in which all energy transitions are possible. The QWIP structure 
in this thesis was designed to detect a laser wavelength of 1.06 |im and a wavelength in 
the 8-10 |im atmospheric window. 

The relevance of a detector that is tuned to these wavelengths is that it can be used 
on military aircraft as a laser spot tracker and an infrared imager providing much greater 
accuracy and dependability than older systems. 
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I. INTRODUCTION 


A. QUANTUM WELL INFRARED PHOTODETECTOR (QWIP) 

The past twenty years have seen an explosion in the realm of infrared detection 
technology fueled by improvements in III-V compound semiconductor technology and by 
new semiconductor growth methods such as molecular beam epitaxy (MBE) and metal 
organic chemical vapor deposition (MOCVD). These technologies have made it possible 
to create artificially structured superlattices for use in novel device applications. The 
QWIPs operate by absorbing incoming photons and exciting carriers from the ground 
state to an excited state, which is known as intersubband transitions and is shown 
schematically in Figure 1.1. These excited carriers, or photoelectrons, can then be used 
to create a photocurrent through the application of an external bias. 



Figure 1.1 Diagram showing quantum well structure, energy levelsenergy 

gap and photon absorption {fico). 

QWIPs have an advantage over other infrared detecting materials such as 

HgCdTe (MCT), which is shown in Figure 1.2, because they have much larger bandgaps 

and therefore are stronger and easier to manufacture. QWIPs also have an advantage 
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over MCT devices due to the fact that MCT devices utilize interband transitions, which 
limit the ability of the detector to be “tuned” to a specific wavelength. These advantages 
have made QWIPs the focus of recent infrared detection research. 



Figure 1.2 Diagram of HgCdTe bandstructure showing the energy gap and the photon 

absorption across the gap. 

The idea of using semiconductors to create quantum wells was first proposed by 
Esaki and Tsu (1969) and since then many studies have been done on artificial quantum 
wells, which has led to a detailed understanding of energy-level spacing, intersubband 
and interband selection rules, line shapes and intersubband absorption. The earliest 
working QWIP was designed using GaAs/AlGaAs and was based on bound-to-bound 
transitions (Levine et al. 1987). The responsivity of this first device was extremely low 
because the photoelectrons had to tunnel through large barriers and therefore the device 
had to operate under a CO 2 laser. This low responsivity led to the use of bound-to- 
continuum transitions, which was first proposed by Coon et al. (1984). The first bound- 
to-continuum QWIP was demonstrated to have a much higher responsivity and a larger 
spectral width (Hasnain et al, 1989), but had a small peak absorbance due to the 
continuum of states above the well. 
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Further advances in QWIP design have led to structures that are based on bound- 
to-quasi-continuum transitions (Levine et al, 1991), bound-to-miniband transitions (Yu 
et al, 1991) and bound-to-quasi-bound transitions (Gunapala et al, 1996). These more 
recent configurations have shown good detector performance due to the fact that in all of 
these cases the ground state electrons do not flow in response to an external bias while 
the photoelectrons can create a photocurrent with the use of a relatively small bias. 
Multiple structural optimizations have also been made that improve QWIP performance. 
Increasing the barrier size has reduced the tunneling current by many orders of magnitude 
(Levine et al, 1991), lowering the excited state from the continuum into the quasi-bound 
region has been shown to reduce the dark current from thermionic emission by a factor of 
~12 at 70K (Gunapala et al, 1996) and adding a grating on top of the device to increase 
the electric field polarization normal to the quantum wells in order to increase the 
absorption strength (Karunasiri et al, 1995). All of these QWIPs have been shown to 
work well as infrared detectors, but they are limited to the detection of a narrow band of 
wavelengths. This limitation in QWIP technology has recently been overcome through 
the design of “multicolor”, or multiple-wavelength, QWIPs. 

Multicolor QWIPs are important because the peak wavelength emitted by an 
object shifts as its temperature changes. This shift in wavelength can cause performance 
limitations in single color QWIPs due to the peak radiation wavelength being far 
removed from the peak detection wavelength of the device, which will limit the 
discrimination ability of the detector. One way in which to improve this limitation is to 
design the QWIP to be sensitive to several wavelengths so that the emissivity changes of 
an object can be detected. 

Some of the earliest methods of creating multicolor QWIPs were to build several 

stacks of square quantum wells, each with a different peak detection wavelength (Kock et 

al, 1992, Grave et al, 1992, Liu et al, 1993). The multicolor response of these stacks 

was achieved by attaching a contact to each stack (Kock et al, 1992) or by controlling 

the external bias across the device to activate each stack sequentially (Grave et al, 1992, 

Liu et al, 1993). Another method of creating a multicolor device was to use bandgap 

engineering technology to create an asymmetric quantum well, where the transitions from 

the ground state to several excited states are allowed and occur at differing wavelengths 
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(Berger et al., 1992, Martinet et al, 1992, and Chiang et al, 1996). The main advantage 
of using this system is that it requires only one set of quantum wells to aehieve multieolor 
detection, but the downside is that in order to get the photoelectrons from the lower 
excited states out of the well a large bias has to be applied to the device, which increases 
the tunneling current. One method of overcoming this large tunneling current is to create 
a superlattice in which an asymmetric well is surrounded by multiple square wells so that 
the first miniband is aligned with the lowest excited energy level in the asymmetric well 
(Ting et al, 1997). The advantage of this system is that photoelectrons in the lowest 
excited energy level have a higher mobility due to the presence of the miniband, and 
therefore the responsivity of this structure is greatly improved, but the complexity of the 
design makes fabrication more difficult. 


B, PURPOSE OF THIS THESIS 

The purpose of this thesis is to design an asymmetric quantum well capable of 
detecting a laser wavelength of 1.06 |im and a middle infrared wavelength in the 8-10 |im 
range through the use of both intersubband and interband transitions. The scope of this 
research will cover several aspects of semiconductor technology including theoretical 
study of band structure, intersubband and interband transition selection rules, 
intersubband absorption, structural design of the quantum well, and modeling of the 
effects of an external bias on the quantum well. 

The solution to the Schrodinger equation in a QWIP with finite potential heights 
cannot be solved analytically, but must be calculated using techniques such as the transfer 
matrix method (Vassel et al, 1983, Lui et al, 1986, Wang et al, 1989). The transfer 
matrix method utilizes the continuity and boundary conditions of the wavefunctions to 
obtain the energy eigenvalues and wavefunctions numerically. These values are then 
used to determine the absorption spectrum of a device, which is useful in optimizing the 
device performance. 

In order to achieve two-color detection, it is necessary to design a quantum well 
structure having two energy transitions corresponding to the two detection wavelengths. 

In a square quantum well, this is difficult because intersubband transitions can only occur 
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between states with opposite parity aceording to the seleetion rules. In addition, the 
transition probability drops rapidly as the quantum number of the exeited state inereases. 
Introducing an asymmetry into the quantum structure and thereby altering the parity of 
the wavefunctions and allowing all transitions to occur can overcome this restriction 
(Martinet et al, 1992). The quantum well structure designed in this thesis makes use of 
both intersubband and interband transitions to achieve multicolor detection. This method 
differs from those studied in the past, which utilized only intersubband transitions to 
obtain multicolor detection (Kock et al, 1992, Grave et al, 1992, Liu et al, 1993, Ting 
et al. 1997). The advantages of the method used in this thesis are that the use of 
interband transitions allows for much shorter wavelength detection while the longer 
wavelengths can be detected using intersubband transitions. Furthermore, the asymmetry 
of the quantum well allows the excitation of electrons from the ground state to the first 
excited state in the conduction band. This allows for the extraction of the photoexcited 
electrons to occur at a lower bias, which reduces the dark current dramatically. 


C. MILITARY RELEVANCE 

A multicolor QWIP is relevant in the military today because it can be used as both 
an infrared imager and a laser spot tracker for use in laser-guided weapons delivery. The 
quantum well design utilized in this thesis is capable of detecting wavelengths as low as 1 
|im and could easily be tuned for use with any combat laser designator wavelengths. In 
addition to the 1.06 |im laser wavelength, this quantum well is designed to 
simultaneously detect infrared radiation in the 8-10 |im range, supporting infrared 
imagery. 

The applications of infrared and laser detection technology in the military are 
legion. Specific benefits of QWIP detector devices are immediately evident in the 
following case: Military airborne weapon systems often utilize a FLIR (Forward 
Looking InfraRed) for target acquisition and tracking and a laser target designator/range 
finder that are often co-located in a single weapons pod (e.g., LANTIRN) (Hewish et al, 
2002). Some aircraft, such as the F-18, also employ a laser spot tracker, which is located 
in a separate, externally mounted pod. The addition of the separate laser spot tracker 
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(LSI) provides the pilot with a visual laser designation eue, whieh is espeeially useful 
when designation oeeurs from an offboard souree (e.g., Speeial Ops, “Buddy lasing”, 
ete.). The LST affords much greater accuracy and confidence in ordinance delivery 
while expanding employment options. However, boresight errors are inherent in a multi¬ 
pod system and degradations due to misalignment are common. Another disadvantage is 
the loss of an external weapons station for carriage of a sensor. 

The quantum well detector design proposed in this thesis will provide an 
advantage to both airborne systems, which will lead to better accuracy and lethality. In a 
system such as LANTIRN, it will provide laser spot tracking capability where one 
currently does not exist. For aircraft such as the F-18, the introduction of a laser spot 
tracker into the weapons pod will not only eliminate boresight errors but also free up the 
weapons station often occupied by the laser spot tracker, resulting in lethality and 
survivability benefits. 
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II. QUANTUM WELL STRUCTURE AND INTERSUBBAND 

ABSORPTION 


A. INTRODUCTION 

Intersubband absorption is an important aspect in the design of an infrared 
detector and has been studied extensively using multi-band methods (Shik, 1992, Peng et 
al., 1995, Flatte et al, 1995) and also using band-mixing (Peng et al, 1995). These 
methods for studying absorption are very sophisticated and may not be suitable for use in 
designing QWIPs. When quantum well structures are formed using large bandgap 
materials, the envelope function approximation can be used instead to provide a more 
simple approach to determining design parameters. In the envelope function 
approximation, the Schrodinger equation is solved by using the transfer matrix method, 
which has been discussed by several researchers (Vassel et al, 1983, Lui et al, 1986, 
Wang et al, 1989) and has become the method most commonly used to model quantum 
well technology today. This chapter will cover the transfer matrix method, band 
structure, selection rules for energy transitions and intersubband absorption. 


B, TRANSFER MATRIX METHOD 


The Schrodinger equation inside an infinite square potential well, which is shown 
in Figure 2.1, is given by: 


2m^ dz^ 




( 2 . 1 ) 


where is the mass of the electron, fi is the reduced Planck’s constant, z is the growth 
direction and ij/^ andE’^ are the wavefunction and energy eigenvalues, respectively. 

Finding the solution to (2.1) is relatively easy due to the conditions that the wave 
functions must go to zero at the edges of the well and that the electron must be within the 
well. In this case the quantized energy levels have the values (West et al., 1985): 
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( 2 . 2 ) 


fl^TT^rP' 

^ forn= 1,2,3... 

2m/ 

where L is the width of the well. The wavefunetions in this ease are given by 


[2 . 

f TTnz^ 

— sm 


u 

1 L J 


(2.3) 


where the , |— term is known as the normalization faetor and eomes from the eondition 
L 

that 


j \¥ni^) f dz = \. 


(2.4) 



Figure 2.1 Sehematie drawing of an infinite square quantum well with energy levels. 


In the ease of an infinite potential well the approaeh used above is relatively 
simple, but it beeomes tedious and diffieult to solve when multiple layers are added as in 
the ease of a finite step quantum well (Yuh et al, 1989). In this seetion the transfer 
matrix method will be diseussed and shown to be able to handle potential wells with 
many layers. 

In a heterostrueture, the quantum well is formed by eombining semieonduetor 
materials that have differing bandgaps. The simplest form of this is a square quantum 
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well in which material A is surrounded by material B, which has a higher bandgap and 
which is thick enough not to allow for tunneling of the confined electron. In this case the 
electrons will be confined to the smaller bandgap region and using the envelope function 
approximation (also known as the effective mass approximation) the energy eigenvalues 
can be derived (Kane, 1957). 

The envelope function approximation makes use of two assumptions; (1) that the 
potential difference between material A and B at the interface is localized (i.e. step 
function) and (2) the interface potential does not mix the band-edge wavefunctions, but 
shifts them. With these assumptions the wavefunctions take the form (Coon et al, 1984) 

T(r) = M(r)e''^-'’<"’V(2) (2.5) 


where ^(r) is the bandedge cell periodic function, z is the crystal growth direction, 

^ikp(x,y) -g plane-wave perpendicular to z, p(x,y) is the direction vector in the xy 
plane, k is the transverse wave vector, and ^^(z) is the envelope wave function that 
satisfies the Schrodinger equation 




2m (z) dz‘ 


- + V(z) 


y/{z) = Ey/{z) 


( 2 . 6 ) 


where m*{z) is the effective mass in the direction of growth and V[z) is the potential 
height. 


The transfer matrix method utilizes the piece-wise constant potential V[z) to 

simplify much of the calculations of the envelope function approximation so that the 
wave functions in each layer are given as a linear combination of a right-traveling wave, 
, and a left-traveling wave, , where is the wave vector in the n* region (Note 
that is real in the layer where E^>V (z) and imaginary otherwise). Using this fact 
and applying the boundary conditions 


¥ni^n) = ¥nA^n) 


1 ^¥ni^n) 


1 






(2.7) 
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we can mathematically determine the energy eigenvalues. 


V(z) 



Figure 2.2 Arbitrary potential well of a multi-layered quantum well structure. 


Figure 2.2 shows an arbitrary quantum well structure in which the thickness of 
each layer is given by and m* is the effective mass of the electron in each 

region. The wavefunctions for each region are given by 




-ik„z 


n = 0,l,2,...,n,n+l 


( 2 . 8 ) 


where and are the coefficients of the left and right traveling wave in the n*’’ region 
and the wave vector in the n'*' region is defined as 



(2.9) 


where is the potential height in the n**' region. By applying the boundary conditions a 

relationship can be developed between the coefficients of the wavefunctions in 
neighboring regions. 

n= 1,2,...,N-1 (2.10) 


^n+l 

= M„ 


n 


where is the transfer matrix for the interface at z^ 
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( 2 . 11 ) 


M 


n 


1 (1+X,)e'‘*" 


(l+Z,)e“'^*"“*"+‘^^" 


y *.'"..L 

/C« , * 


( 2 . 12 ) 


From this formula a relationship can further be found between the first and the last region 
in the quantum structure 


A 


= M^_j...M2Mi 



= M 



mil 

m,2 

~A 





_m2i 

m22_ 

A. 


(2.13) 


where M is the composite transfer matrix for the entire structure and is given by 


M = M^_i...M2Mi 


m,i mi2 


(2.14) 


If the potential heights in the two outer regions are higher than those in between them 
there will be bound energy states. For these bound states the wavefunctions must go to 
zero at ±oo and this requirement means that =0 and Bjy = 0 : 


1 

1 


1 

s 

mi2 

1 

o 

1_ 

1 

O 

1 _ 


_J 

-1 

<N 

1- 

> 

1_ 


which is satisfied when 


(2.15) 


m22(B) = 0 (2.16) 

Equation (2.16) is then tantamount to an eigenvalue condition on the allowed energies. 
We wrote a program to plot m 22 {E) versus E and locate the zeroes of the function, which 
determine the energy eigenvalues. The zeroes of m 22 {E) are found by searching for 
values of E such that m 22 {E)\m 22 {E + EE)\ < 0 ; in this way we know that the root lies 

between E and E + EE . The accuracy of the eigenvalue is determined by EE ; for most 
cases a value of EE = 1 meV was used. 
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An example of the transfer matrix method for a square quantum well is shown in 
the next seetion, where the energy levels for the quantum well shown in Figure 2.3 are 
derived. 

V(z) 


m 


V. 


m 


V2 = 0 


V3=V, 


m*3=m*^ 


kg k.| 


z 


0 


Figure 2.3 Schematic diagram of a finite square well with parameters used for 

calculation. 


The composite transfer matrix M for this quantum well can be written as 


M = M2M,=i 




l+Zi l-Z 
I-Zl I+Zl 


(2.17) 


2m 2m 

where k^=k^= iK = i. —~ k 2 = ^ , so that, 

\ h \ h 


m 


22 ' 


i{k2+k2)a 




^-i{k2-k2)a 


(2.18) 


The bound states occur when m22iE)=0. Figure 2.4 shows a plot of ^ function 

of energy. 
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m22(E) vs. Energy 



Figure 2.4 Plot of versus E for InP !well with a = 150 A. 

Although only the conduction band was considered in this case, the valence band follows 
the same procedure and the energy eigenvalues in the valence band can be solved for 
with the same program keeping in mind that the effective masses and potential heights 
will differ in the valence band. Figure 2.5 shows the conduction band energy levels and 
wavefunctions of the quantum well in Figure 2.4. 
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Conduction Band Energy Levels and Wavefunctions 



Figure 2.5 Conduction band energy levels and wavefunctions of the quantum well in 

Figure 2.4. 

C. SELECTION RULES FOR OPTICAL TRANSITIONS 

In a quantum well structure made of semiconductors, the wavefunctions consist of 
envelope and Bloch functions. This means that transitions between both Bloch and 
envelope states are allowed within the structure. Transitions between Bloch states 
correspond to interband transitions, while transitions between envelope states correspond 
to intersubband transitions. The transition rate for either interband or intersubband 
transitions is given as (Coon et al, 1984) 

W = — Xl< I I > 1 ' S[E^-E,-no)\ (2.19) 

n ij 

where E. andE'^ are the initial and final energy states and hO) is the incoming photon 

energy. The initial and final wavefunctions in the transition rate are of the form of 
Equation (2.5), 
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( 2 . 20 ) 




The interaction potential is defined in the dipole approximation as (Coon et al, 1984) 



in 


J/2 


Is^n^CiX: 


A. 

e p 


( 2 . 21 ) 


where / is the incident photon flux, c is the speed of light, e is the electron charge, £q is 
the electric permittivity, is the index of refraction, a> is the angular frequency of the 

incoming photon, e is the unit vector of the electric polarization of the incident light and 
p is the momentum operator. The integral in the transition rate is known as the transition 
matrix element, which is defined as 

=< M^(r)e I Vp I >< y/f{z) \ > (2.22) 

+ < I Vp I X M^(r) I M,(r) >. 


For the case of interband transitions ( if/^ = y/^, Wf-¥c) second term in the transition 

matrix element is zero due to the orthogonality of the Bloch functions in the conduction 
and valence bands. On the other hand, the first term is zero for intersubband transitions 
(\|/^ =y/y. If/j = 11 / 2 ) due to the fact that the Hamiltonian used to obtain xj/p and if/^ is 

Hermitian. Taking this into account, the transition matrix elements for both transitions 
are as follows 


M,, =< M,(r) I I M.(r) X y/^{z) \ Xf/^V) > 

Mi2 =< 1^2(^) I Vp I Xf/^{z) >< u^ir) I Mi(r) > 


(2.23) 


where the term, giving a delta function for the conservation of the xy plane 

momentum, has been dropped. It is easy to see that the only way to obtain non-zero 
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transition matrix elements is for both < y/^{z) \ 0 and < y/y{z) \ \ 0 • 

For a square quantum well this implies that only even-to-even or odd-to-odd interband 
transitions are possible, while for intersubband transitions, and only even-to- 

odd or odd-to-even transitions are possible as shown in figure 2.6. 



Figure 2.6 Sehematic diagram of a square quantum well showing allowed and 

forbidden (X) transitions. 

If a step is added within the well in Figure 2.6 then the wavefunctions are no 
longer only even or odd, but a combination of the two. In this case the transition matrix 
elements will be non-zero for all transitions as is shown in figure 2.7. 
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g 


Figure 2.7 


t 



Schematic diagram of a step quantum well showing that all transitions are 

allowed. 


This step is vital to this thesis because it allows for the excitation of electrons from 
the and states to the E^ state to be collected as photocurrent. The optimization of 
transition wavelengths and the transition probability are discussed in the next section. 


D. INTERSUBBAND ABSORPTION 

Absorption is an important factor in the design of a multicolor infrared detector 
because it is one of the only measurable quantities of a sample. The absorption 
coefficient is proportional to the amount of incident light that is absorbed in the sample 
and therefore should be high in order to achieve an efficient detector. This thesis will 
include absorption from both bound-to-bound and bound-to-continuum transitions, which 
are shown in figure 2.8. 
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(a) (b) 

Figure 2.8 Diagram of (a) bound-to-continuum and (b) bound-to-bound transitions. 


For the case of bound-to-bound transitions there is only one initial state and one 
final state, but there are N = N^LA electrons occupied in the initial state due to freedom 
of motion of the electrons in the xy plane, where N is the number of electrons in the 
initial state, is the doping concentration in electrons per volume, L is the length of 

the well and A is the area of the well. Since all of these electrons have the same energy 
along the growth direction the sum in Equation (2.19) is over just one initial and final 
state multiplied by N . The delta function in Equation (2.19) can be approximated due to 
finite lifetime of the excited state by the normalized Eorentzian function given as (Yariv, 
1997, pp. 162) 




(2.24) 


where E is the full width half maximum value of the electron in the state . Elsing the 
interaction potential from Equation (2.21) and the momentum operator 




(2.25) 


the transition rate can be found as 


l< r, 14-1 ¥, >t g(*®) 


(2.26) 
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where < | | 'Pj >~< | | > for intersubband transitions and 


C 


B-B 


N,e^n^ 

[m) e^n^c 


(2.27) 


For the ease of bound-to-eontinuum transitions there are still N initial states with 
energy E ., but now there is also a eontinuum of final energy states and the summation 
over these final states amounts to replaeing the delta funetion with the joint density of 
states p{E^), whieh is given by (Choi, 1997, pp. 131-133) 




27r^n\E.-V,) 


(2.28) 


where ml is the barrier effeetive mass and FJ, is the barrier potential height. Therefore 
the transition rate for bound-to-continuum states is defined as 


W, 


B-C 




2m, 


inojfiE.-V,) 


\< ¥f 


A. 

dz 


¥i >r 


(2.29) 


where 


2+,2 




I =*= \2 

2(m j epi^c 


(2.30) 


The absorption eoefficient a{hco) is defined as 


a{hQ)) = 


number of transitions per unit volume and time 
ineident photon flux 



(2.31) 


and the absorption eoefficient for both bound-to-bound and bound-to-continuum 
transitions can be derived. 


CCB-Bi^co) = \<¥f\^\ ¥i >f gi^O)) 

nO) dz 


2m 

a„ AfiO)) = C„ r I -^— - \<¥f\ — I ¥i >1 


d I i2 

— I w, >1 

dz 


(2.32) 
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An example of the absorption coefficient for bound-to-bound and bound-to-continuum 
transitions is shown in figures 2.9 and 2.10. The sharp peak in figure 2.9 results from the 
bound-to-bound transitions occurring whenever the incident photon has energy equal to 
E^-E^, whereas the shape of the absorption curve in figure 2.10 is relatively broad due 

to the continuum of energies above the well. The amplitude of the bound-to-continuum 
absorption coefficient decreases as the photon energy increases because the density of 
states decreases the further above the well you go. 


Bound-to-Bound Absorption 



photon energy (eV) 


Figure 2.9 Bound-to-bound absorption curve for (Ino. 53 Gao. 47 As)o.i 5 (InP)o. 85 / 
Ino.53Gao.47As /(Ino.53Gao.47As)o.45(InP)o.55 with a=25 A and b =44 A with no external bias 

and r = 10 meV . 
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Bound-to-Continuum Absorption 



Figure 2.10 Bound-to-continuum absorption curve for the quantum well in Figure 2.9 

with no external bias. 


E. SUMMARY 

In this chapter the transfer matrix method was developed and used to calculate the 
energy levels and wavefunctions in a quantum well. Selection rules for interband and 
intersubband transitions and their importance in detector design were discussed and the 
absorption coefficient for both bound-to-bound and bound-to-continuum transitions was 
derived theoretically using Fermi’s golden rule (2.19). In the next chapter, we will 
discuss some of the design aspects that are considered when designing a multicolor 
QWIP and how the addition of an external bias will effect these considerations. 
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III. DESIGN OF A QW STRUCTURE FOR INFRARED 

DETECTION 


A. INTRODUCTION 


As discussed in Chapter II, surrounding a low-bandgap material with a high- 
bandgap material will form a quantum well. The choiee of which materials to use is 
important as it will affect the transition wavelengths as well as the absorption coefficient, 
and therefore the effieiency of the detector. Some of the factors to consider in the ehoice 
of materials to use are bandgap, lattiee constant, and what wavelengths are to be detected. 
During this thesis multiple materials were tried in order to optimize the structure so that it 
would deteet a wavelength of 1.06 pm and a wavelength of 8-10 pm. No matter whieh 
materials were used in the structure, their bandgaps and effective masses were 
determined by Vegard’s Law with a bowing parameter to take into aecount non¬ 
uniformities of alloying (Singh, 1993, pp. 185). 




= xE^ + (1 - x)E^ + x(l - x)C 


AB 


m. 


^ X 1-x^ 

ml j 


(3.1) 


where x is the molar pereentage of material A and is the experimentally determined 
bowing parameter. An example of how the composite energy gaps, at the L-point, and 
relative effective masses ehange with molar pereentage is shown in Figures 3.1 and 3.2 
where the energy gaps and effective masses for eaeh material are given in Table 3.1 
(Note that the bowing parameter is not considered in these eases). 
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Energy Gap (eV) 


Table 3. 


Material 

Energy Gap (eV) 

Relative effective mass [m* 

AlAs 

2.671 

0.15 

GaAs 

1.424 

0.067 

In As 

0.36 

0.028 


Energy gap and relative effeetive mass values for speeific materials. From; 


Singh, 1993, pp.184-185 
Energy Gap vs. Molar Percentage 



Figure 3.1 Diagram of Energy Gap vs. Molar percentage for 
Alfia^_^As and Ga. Jn^_^As at the T-point. 
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Relative Effective Mass vs. Molar Percentage 



Figure 3.2 Diagram of Relative effeetive mass vs. Molar percentage for 

Alfia^_^As and GaJn^_^As 


The potential heights in the conduction and valence band depend on the difference 
in energy gaps between two regions and are not the same, as might be expected, but differ 
so that the conduction band potential height is larger than the valence band. This 
difference is difficult to predict theoretically and can only be determined for each 
material experimentally, however, when modeling a quantum well structure it is assumed 
that to a good approximation the potential heights in the conduction and valence bands is 
given as 



In early QWIP research, studies focused on quantum well structures that utilized 
bound-to-bound transitions (Levine et al. 1987). In a bound-to-bound QWIP the 
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electrons in the lowest state are excited by incoming photons into the excited state and 
then tunnel through the barrier under an external bias to form a photocurrent as is shown 
in figure 3.3. 



Figure 3.3 Schematic Diagram showing (a) bound-to-continuum and (b) bound-to- 

bound electron transport. 

The barrier thickness was designed to be large enough to prevent tunneling from the 
lowest state (typically about 500 A); however, this barrier thickness reduced the 
photoresponse due to the reduction in tunneling probability of the photoelectrons. 
Another problem with these early designs was that the photoresponse spectrum was 
relatively narrow due to there being only two states and was not practical for infrared 
imaging. These problems ushered in the second generation of QWIPs, which used 
bound-to-continuum transitions (Hasnain et al, 1989). In these structures, electrons in 
the ground state are excited into the continuum where they contribute to the photocurrent 
without tunneling (Figure 3.3) and therefore the barrier thickness has less bearing on the 
photocurrent. Due to the continuum of states and the multiple transitions that can occur, 
the photoresponse spectrum was much wider than in the bound-to-bound design, but the 
absorption coefficient was much smaller due to the low density of states above the well. 
It is easy to see that there is a trade-off between having good transport of photoelectrons 
and strong absorption and that both are important in the design of a QWIP. In the 
following sections, these factors will be used to motivate the design of a quantum well 
structure for two-color infrared detection. 
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B, STEP QUANTUM WELL STRUCTURE FOR DUAL COLOR 

DETECTION 

In the past, attempts have been made to deteet multiple “eolors”, or wavelengths, 
with an infrared imager by using square quantum wells (Koek et al, 1992, Grave et al, 
1992, Liu et al, 1993), but as was shown in seetion 2.3, transition seleetion rules prohibit 
the transition from the first valenee level to the seeond eonduetion level and therefore that 
wavelength eannot be deteeted. This implies that only one wavelength of light ean be 
detected with a square quantum well. One method used to detect more than one 
wavelength of light using square quantum wells was to stack multiple wells on top of 
each other and hold each well responsible for detecting one specific wavelength of light 
(Koek et al, 1992, Grave et al, 1992, Liu et al, 1993). This type of structure is easy to 
design, but very difficult to fabricate due to the need for multiple contacts at each well. 
Recently researchers have been exploring the use of asymmetric, or step, quantum wells 
(Ting et al, 1997) to detect multiple wavelengths. 

In a step quantum well, the symmetry of the wavefunctions is broken and 
transitions that were forbidden in a square quantum well can occur. For this thesis, the 
optimized design of the step quantum well is to have two energy levels in the conduction 
band, one above and below the step, which account for the middle infrared wavelength, 
and one energy level in the valence band, below the step, that accounts for the 1.06 pm 
transition as illustrated in figure 3.4. The second energy level in the conduction band is 
designed to be as near the top of the quantum well as possible in order to get a good 
photocurrent with very little bias. 
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V(z) 



The application of the bias will cause the quantum well to “slant” and will force 
the second energy level into the continuum allowing the photoexcited electrons to be 
extracted without tunneling through the barrier. The energy levels will be affected by the 
bias and therefore the peak detection wavelengths will shift. The asymmetry of the 

cFb 

quantum well will give rise to a linear Stark shift which is approximately given by 

where e is the electron charge, F is the applied electric field in volts per meter and b is 
the length of the step in the quantum well. The accurate shift in detection wavelength is 
important in the design since it is necessary to apply an external bias to a device in order 
to extract the photocurrent. Due to the necessity of using an external bias it becomes 
important to be able to model how the peak detection wavelengths will change in order to 
optimize the QWIP structure, in particular for the 1.06 |im wavelength. The next section 
will introduce a transfer matrix method for a quantum well under a bias and this method 
will be used to model the change in transition wavelengths due to an external bias. 


28 




C. TRANSFER MATRIX METHOD UNDER AN APPLIED BIAS 

When a quantum well structure is subjected to an applied bias the potential well 
will tilt in the direction opposite the electric field and therefore the potential of an 
electron in the quantum well is given as 

V{z) = V„-eFz (3.3) 


where is the potential height of the well, eis the electron charge, F is the 

applied electric field in volts per meter and z is the direction of growth. An example of 
the effect of an external bias on a square quantum well structure is shown in Figure 3.5 



Figure 3.5 Diagram of a square quantum well under an external bias. 
Substituting Equation (3.3) into the Schrodinger equation 


^ 2m* dz^ 


+ V(z) 




(3.4) 


the following formula can be derived 


dz' 


n 


where 


(3.5) 


29 



(3.6) 


Vn = 


E-Vn 

e 


and where is the potential height in the n* region. Equation (3.5) ean further be 
simplified using the following eoordinate transformation 

= (3.7) 


and the result is of the form of the Airy differential equation 






(3.8) 


with 




f;j2^2V3 




(3.9) 


The solution to Equation (3.8) is a linear eombination of Airy functions given as 
(Vatannia et al, 1996) 

(z) = AJi{p^ (z)) + B„Bi{Pn (^)) (3.10) 


where Ai and Bi are Airy functions of the first and second kind, respectively. By 
substituting Equation (3.10) into the boundary conditions given by Equation (2.7) a 
relationship can be developed between the coefficients of the wavefunctions in each 
region. 



= M„ 

A+\ 

Bn'_ 


_Bn+\ _ 


for n = 1,2,3,...A-1 


(3.11) 


Eurthermore, a relationship can be developed between the coefficients in the first and last 
region. 


A 





= M 


_ 

mu 

m,2 


Bn} 


Bn} 



^22 _ 

Bn 


(3.12) 
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where 


M„ =;r 


Aiia„JBi'ia„)- Biia„JBi'ia„)- CT„Bi{a„)BiXa„j' 

yanAi(an)AiXanJ-AiXaJAi{a„J aJi{aJBiXa„J- AiXaJBi{a„J 


(3.13) 


%+i^y 


and 






C(„ 


V y 


'iml.eF''''’ 


i-^n-Vn)=rn(-^n-Vn) 


(3.14) 


/ * \ 
^nYuFl 


(3.15) 


The primes in Equation (3.13) represent the derivatives of the Airy functions and the 

factor of 7t comes from the fact that [Ai{aJBiX(Xj- AiXo(JBi(aj) = —. In the case of 

n 

bound states, the wavefunctions must decay to zero as z ^ ±oo . This constraint forces 
= 0 since the second order Airy function, Bi(aJ , tends to oo as z —> -oo . Similarly, 

Aj^ =0 since Ai{a^) tends to 0 ° as z —> 0 °. These conditions cause Equation (3.12) to 
take the form 


1 

1 


1 

_1 

1 

0 

1_ 

0 

1 _ 


<N 

_J 

m22_ 

1 

_ 1 


(3.16) 


and, just like the derivation of the unbiased transfer matrix method, the requirement for 
eigen states is 

m^2iE) = 0. (3.17) 


We wrote a program to find and plot m 22 {E) versus E and find the energy levels when 

Equation (3.17) is satisfied. Eigure 3.6 shows m 22 {E) as a function of energy and figure 

3.7 shows the calculated energy levels and wavefunctions for the conduction band of a 
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step quantum well under an external bias. The parameters of the quantum well were 
taken from a previously published work (Karunasiri et ai, 1990) so that the results found 
using this program eould be eompared with experimental data. 



Energy (eV) 

Figure 3.6 Plot of m 22 (E) vs. E for GoAsIA iQ -iiGaQ^ 2 As well with a=60 A 

and b=90 A and an external electrie field, = 3.5F///m . 
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Conduction Band Energy Levels and Wavefunctions 

T-1-1-1-1-1-1-r 



I _I_I_I_I_I_I_I_I_ I 

-150 -100 -50 0 50 100 150 200 250 300 

z (angstroms) 


Figure 3.7 The first three conduetion band energy levels and wavefunctions of the 

quantum well in figure 3.6. 

The eigen states and wavefunctions of the step well formed in the valence band 
can be obtained using the same approach as above with corresponding potential heights 
and effective masses. Figure 3.8 shows the calculated energy levels and wavefunctions 
in the valence band for heavy holes. 
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Valance Band Energy Levels and Wavefunctions 



Figure 3.8 Valence band energies and wavefunctions for the step well in Figure 3.6. 

The introduction of the external bias will cause a shift in the energy levels of the 
asymmetric quantum well and therefore, the absorption coefficient defined by Equation 
(2.32) will also change. For the QWIP design used in this thesis, the addition of the bias 
will force the excited state, which was very near the top of the well, into the continuum 
so that only the bound-to-continuum absorption coefficient is considered when a bias is 
introduced. Figure 3.9 demonstrates how the external bias will affect the peak detection 
wavelength. 
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Bound-to-Continuum Absorption 



wavelength (|j.m) 

Figure 3.9 Plot of the bound-to-eontinuum absorption vs. wavelength for the 
(^'* 0 . 53 ^^ 0 . 47 ^'^)o.l 5 (^'^P)o .85 ^ ^” 0 . 53 ^'^ 0 . 47 "^'^ ^(^”o. 53 ^'^ 0 . 47 ^'^)o .55 (^'^P)o .45 Well wlth a — 25 A and b — 
44 A and with a bias of 0 V/|im and 1.5 V/|im. 

This section has introduced a transfer matrix method for a quantum well structure 
under an external bias. This method was used to calculate the energy levels, 
wavefunctions and absorption coefficient for a quantum well, which is useful in 
optimizing detector design. The next section will compare the output of the computer 
model designed in this thesis with experimental data in order to tell the accuracy of the 
computer model. 


D. COMPARISON WITH EXPERIMENTAL DATA 

The quantum well structure that the computer model was compared to came from 
a paper written by Karunasiri et al. (1990). This paper utilized a 
/ GaAs I AIq ^^GaQ^ 2 ^s step quantum well with a = 60 A and b = 90 A. This 

structure has four energy levels in its conduction band, but the paper focused on the 
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bound-to-bound absorption from the 1-2 and 1-3 transitions due to the minute absorption 
strength of the 1-4 transition. Figure 3.10 shows the energy levels of the conduction band 
and Table 3.2 shows the effective masses and energy gaps used in the computer model. 


Energy Levels of the Asymmetric Quantum Well 
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Figure 3.10 Energy levels of the ! GaAs ! well with a = 60A 

and b = 90A. 


Layer 

m* (kg) 

E, (eV) 

Alo,44Gao,56As 

0.1082 

1.9727 

GaAs 

0.067 

1.424 

Alo,i8Gao,82As 

0.0794 

1.6485 


Table 3.2 Parameters used in the quantum well in Figure 3.10 

Since this is a step quantum well all energy transitions are allowed and only the bound- 
to-bound absorption is considered in this analysis to compare with the experimental data. 
Figure 3.11 shows the calculated absorption coefficient for the 1-2 and 2-3 energy 
transitions under both positive and negative external bias. 
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Bound-to-Bound Absorption for 1-2 and 1-3 Transitions 



Figure 3.11 Bound-to-bound absorption curves for 1-2 and 1-3 transitions of the 
quantum well in Figure 4.6 with r=15meF . 

Since the experimental data is measured in absorbance, the absorption coefficient is 
converted to absorbance using the following relation 

Absorbance = x{\oge)a{ho)) = W (\l°)){\oge)a{hcd) 

(3.18) 

= (2.56*10"' cm)(log e)a{n a) 

where W is the width of the absorbing material and the sin^(17°) term comes from the 

polarization of the light hitting the device at the Brewster angle {6 = 13°). Using this 
equation the absorbance can be found for both transitions and is shown in Figure 3.12. 
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Bound-to-Bound Absorbance for 1-2 and 1-3 Transitions 



Figure 3.12 Bound-to-bound absorbance curves for 1-2 and 1-3 transitions of the 
quantum well in Figure 4.6 with r=15meF . 

This Figure can be compared with the absorbance curves obtained experimentally, which 
is shown in Figure 3.13 (Note that the absorbance data in the Figure is plotted in 
decreasing order of wavenumber). 
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Figure 3.13 Bound-to-bound absorbance curves obtained experimentally using the 
quantum well in Figure 4.6. From: Karunasiri et a/., 1990. 


Comparing Figures 3.12 and 3.13 it can easily be seen that the values for the absorbance 
obtained from our model and those obtained experimentally are very close and that the 
computer model used in this thesis is a good approximation to the real world. 


E. SUMMARY 

This chapter introduced some of the parameters that affect the design of a QWIP 
and a method for detecting multiple infrared wavelengths using an asymmetric quantum 
well design. The optimum design characteristics used in this thesis were also presented 
along with a motivation for modeling the effects of an external bias on a quantum well 
structure. A transfer matrix method utilizing Airy functions was derived to model the 
effects of an external bias and this method was used to show how the peak detection 
wavelength shifts due to the bias. The data obtained using our computer model was then 
compared to experimental data and shown to be a good approximation to real life, 
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granting greater confidence in any information obtained using the model. The next 
chapter introduces the optimum design of the QWIP structure used in this thesis along 
with an analysis of how it will react to an external bias. 
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IV. DESIGN PARAMETERS OF THE TWO-COLOR DETECTOR 


A, INTRODUCTION 

As was mentioned in ehapter III, many different material systems were eonsidered 
in designing the asymmetrie well used in this thesis. At first, a 
Alfia^_^As / GaAs / AlyGa^_yAs step well was eonsidered due to the abundanee of 

experimental data on these two materials and because they are lattice matched, which 
would make them much easier to grow. This structure could be designed to detect the 
middle infrared wavelength (MIRW), but could not detect the 1.06 |im wavelength due to 
its large bandgap. This prompted the use of materials with smaller bandgaps, which 
could be identified by using Figure 4.1 (Tu et al, 1992, pp. 158). 

Lattice mismatch, f 



0.55 0.60 0.65 

Lattice constant (nm) 

Figure 4.1 Diagram of bandgap vs. Lattice constant for different materials. 
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The next GaAs-hz.sQd quantum well design that was tried eonsisted of 
GaAs / In^Ga^_^As / InyGa^_yAs and eould be designed to deteet both the MIRW and 1.06 

|im wavelengths, but not at the same time. The reason that both wavelengths eould not 
be deteeted at the same time using these materials was that in order to deteet the 1.06 |im 
wavelength the energy gaps had to be widened, but this eaused the two eonduetion band 
energy levels to beeome closely spaced, which pushed the MIRW outside of the 8-10 |im 
design window. Another problem with the design of this quantum well was that if the 
molar percentage of Indium became large there was a significant lattice mismatch 
between the two materials, which would introduce dislocations into the sample. In 
order to lessen the lattice mismatch, a quantum well consisting of 
InP ! In^Ga^_^As / InyGa^_yAs was designed because the energy gap of InP and GaAs are 

about the same and because the two materials have much closer lattice constants to 
provide for easier growth. This quantum well structure could also be designed to detect 
the MIRW and 1.06 |im wavelengths, but the problem remained that both wavelengths 
could not be detected at the same time because the energy gap of InP was still too large. 
Following Figure 4.1 this left only one choice of possible materials for the final design: 
the quaternary InGaAsP and As . These materials provided the advantage of 

being lattice-matched and also by varying the amount of InP in the quaternary the 
bandgap could be brought down to the point where both wavelengths could be detected. 
The following section will discuss the quantum well that was designed using these 
materials and will show the theoretical results for the absorption with and without an 
external bias. 


B, OPTIMUM QUANTUM WELL DESIGN 

One of the problems with using the quaternary material in the design of the 
quantum well was that very little experimental data was readily available and no 
theoretical data could be found. However, a paper written by Vurgaftman et al. (2001) 
provided all of the necessary experimental data not only for the quaternary material, but 
also for every other material that was used throughout this thesis. Using the values from 
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this paper, an optimized quantum well was designed with the following layer structure: 
(/Hq 53 ^ 6 ( 047 ^ 45)0 j 5 (/n/’)og 5 //«(, 53^00 47 ^ 45 /(/no 53 Gao 47 ^ 45)0 55 (/nP),, 45 with a thickness of 
the well of 25 A and of the step of 44 A. Figure 4.2 shows the optimized quantum well 
structure, energy levels, and transition wavelengths and Table 4.1 shows the values of the 
effective masses and energy gaps in each layer of the optimum design. 


Energy Levels of the Optimized Quantum Well 
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Figure 4.2 Optimized quantum well structure and energy levels with no external bias. 


Layer 

m (kg) 

Kh (kg) 

E, (eV) 

(lno.53Gao.47As)o.i5(lnP)o.85 

0.0676 

0.4905 

1.3159 

lno. 53 Gao. 47 As 

0.0365 

0.3403 

0.8161 

(lno.53Gao.47As)o.55(lnP)o.45 

0.0519 njg 

0.4242 

1.1181 


Table 4.1 Parameters used to design the final quantum well structure used in this 
thesis, mlf^ is the heavy-hole effective mass and is the electron mass. 

It is important to mention that the values for the energy gaps found in Table 4.1 are the 

values when the temperature, T = OK . The reason for choosing these values over the 
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room temperature {T = 290K ) values is that in order to limit the dark current of the 
device under operating conditions it must be cooled down io T < 11K . The difference in 
the energy gaps and the energy levels of this structure at T = OK and T = 11K is small 
so that for design purposes it was assumed that the device is operating at T = OK. 
Figures 4.3 and 4.4 show the effects of external bias on the quantum well structure and 
energy levels. 


Energy Levels of the Optimized Quantum Well 



Figure 4.3 Optimized quantum well structure and energy levels with Fi,y^t=\.5Vl^im . 
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Energy Levels of the Optimized Quantum Well 



It is easy to see in these Figures that the seeond energy level has indeed been pushed up 
into the continuum where it can easily create a photocurrent. 

The main factor to consider in the evaluation of this design is how the application 
of the external bias will affect the quantum well’s ability to detect the designated 
wavelengths of 1.06 |im and 8-10 |im band. Figure 4.5 shows a comparison between the 
absorption coefficient for the unbiased well where bound-to-bound absorption occurs and 
for the case when =\.5Vl/ilm where bound-to-continuum absorption occurs. 
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Bound-to-Bound and Bound-to-Continuum Absorption vs. Wavelength 



Figure 4.5 Comparison of the absorption eoeffieients for the unbiased well (bound- 
to-bound) and the biased well (bound-to-continuum) with Ff,y^i=\5Vj^im and r=10ff!eF . 

This Figure demonstrates the relative strengths of the bound-to-bound and bound-to- 
eontinuum absorption eoeffieients. The bound-to-bound absorption eoeffieient has a 
narrow speetral width eompared to the broad bound-to-eontinuum speetrum as illustrated 
in Figure 4.6 using the normalized absorption eoeffieients. 
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Normalized Absorption vs. Wavelength 



wavelength (|j.m) 


Figure 4.6 Comparison of the normalized absorption eoeffieients for the unbiased 
well (bound-to-bound) and the biased well (bound-to-eontinuum) with Fg^i=\.Wjiim and 

r=10meV . 

Although the shift from bound-to-bound to bound-to-continuum absorption has dropped 
the relative height of the absorption spectrum it has not dropped to a point where it would 
be difficult to measure. The broadening of the absorption spectrum has also ensured that 
the quantum well structure will still be able to detect the MIRW reasonably well and as 
long as the external bias is not too large, the QWIP should operate as designed. It is 
important to mention that this thesis has not considered interband absorption due to the 
difficulties in defining and integrating the Bloch functions in each material, but it can 
safely be assumed that the interband absorption curve will follow the same pattern that 
the intersubband absorption curve did and therefore the quantum well structure will also 
detect the 1.06 |im wavelength reasonably well. 
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C. SUMMARY 

This chapter introduced the optimized design of the quantum well structure used 
in this thesis along with some of the reasons why the materials used in this quantum well 
were chosen over other materials. It was also shown that the optimized quantum well 
structure would operate as a multicolor detector for the specified wavelengths with or 
without an external bias. 


48 



V. CONCLUSION 


An asymmetric quantum well device capable of achieving two-color detection 
through the use of bandgap engineering was demonstrated and the properties of this 
device were derived and discussed. The ability of the optimized device to detect both 
short infrared wavelengths and middle infrared wavelengths could be useful as a dual 
infrared imager/laser spot tracker in military aircraft for the use of accurate ordinance 
delivery. 

The detailed quantum well structure and eigen energy levels of the device were 
derived through the use of a transfer matrix method with effective mass approximation 
and the intersubband absorption spectrum for both bound-to-bound and bound-to- 
continuum transitions was derived and utilized to optimize the device parameters. A 
transfer matrix method using Airy functions was derived to model the effects of an 
external bias on the device and the output of this computer model was compared to 
experimental data and shown to be a good approximation to real life. The reasoning for 
the choice of the materials used in the optimized design was discussed including a 
discussion of why other materials that were tried did not work for the chosen 
wavelengths. 

This thesis has shown that an efficient laser spot tracker/infrared imager quantum 
well device is plausible and further study is warranted. A sample of the optimized design 
should be grown and tested in order to compare experimental results with the theory 
presented in this thesis. 
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APPENDIX. [PROGRAMMING NOTES AND SELECTED 

PROGRAMS] 


Since many of the calculations presented in this thesis involved a substantial 
amount of matrix operations the choice of Matlab to ereate the computer model was 
manifest. Matlab proved to be an extremely efficient and reliable tool for deriving the 
energy levels, wavefunctions and absorption spectra for the unbiased quantum well 
devices. However, the use of Matlab to model the biased quantum well had some 
limitations due to the Airy funetions. Matlab contains the numerieal values of the Airy 
functions, but as expected, these values become unmanageable when the argument in the 
Airy functions becomes very large or small. An attempt was made to manually insert the 
asymptotic versions of the Airy funetions into the program to overcome this limitation, 
but this too proved to be unmanageable. Therefore the program written to model the 
quantum well under an applied field is limited in its scope. The following is a selection 
of programs used in this thesis to obtain the theoretical data discussed above. 
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%This program is titled aqw_InGaAsP% 
%Created: January 2002 
%Author: Kevin Lantz 


%This is a program to find the energy levels in a two-step quantum well. 

%There are three interfaces in this well and therefore will be three transfer matrices to calculate. 
%The basic design of the well is as follows: 

%Region 1 (-Inf->0)->InGaAsP 

%Region 2 (0->a)—>InGaAs 
%Region 3 (a->L)->InGaAsP 
%Region 4 (L->Inf)—>InGaAsP 


clc 

global E 

% inputs- 

XI = input('Percentage of InGaAs (0-1) in region 1 '); 

X2 = input('Percentage of In (0-1) in region 2 '); 

X3 = input('Percentage of InGaAs (O-I) in region 3 '); 

widthi = Inf; %Extends from -Infinity to 0; 
width2 = input('Length of 0->a in angstroms: '); 

width3 = input('Length of a->L in angstroms: '); 

width4 = Inf; %Extends from L to Infinity; 


%Energy gap values in each region- 

egl = 1.4236*(1-X1) -t 0.81612*X1 - 0.13*(X1*(1-X1)); 
eg2 = 1.519*(1-X2) -t 0.417*X2 - 0.477*(X2*(1-X2)); 
eg3 = 1.4236*(1-X3) -t 0.81612*X3 - 0.13*(X3*(1-X3)); 
eg4 = 1.4236*(1-X1) -t 0.81612*X1 - 0.13*(X1*(1-X1)); 

%Potential Heights in Conduction Band- 

VI =0.6*(egl - eg2); 

V2 = 0; %Choose the origin as the bottom left comer of region 2. 

V3 = 0.6*(eg3 - eg2); 

V4 = V1; 

% - 

y = 1; %Initial value of y to correct the logic of the later condition; 
a=l; 

%constants- 

hbar = 1.055e-34; %hbar in ev*m 

mO = 9.1 le-31; %Rest mass of the electron in eV 

q= 1.602e-19; 

L1 = width 1 * 1 e-10; %Change input of angstroms->meters 

L2 = width2*le-10; 

L3 = (width2-t-width3)*le-10; 

L4 = width4*le-10; 

%Effective electron mass in each region- 

m_effl = l/((Xl/(0.0365*m0)) ((l-Xl)/(0.0795*m0))); 

m_eff2 = l/((X2/(0.026*m0)) ((l-X2)/(0.067*m0))); 

m_efG = l/((X3/(0.0365*m0)) ((l-X3)/(0.0795*m0))); 
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m_eff4 = l/((Xl/(0.0365*m0)) + ((l-Xl)/(0.0795*m0))); 
for E = 0.0001:0.0001:(V1-0.0001); 


%Determine the values of k- 

kl = sqrt((2*m_effl*(E - Vl)*q)/(hbar^2)); 
k2 = sqrt((2*m_eff2*(E - V2)*q)/(hbar^2)); 
k3 = sqrt((2*m_eff3*(E - V3)*q)/(hbar^2)); 
k4 = sqrt((2*m_eff4*(E - V4)*q)/(hbar^2)); 

%Determine the transfer matrices- 


xl = 0; %Interface Locations; 
x2 = L2; 
x3 = L3; 

mil =((1/2)*((1 +((kl*m_eff2)/(k2*m_effl)))*exp(i*(xl*(kl - k2))))); 
ml2 = ((1/2)*((1 - ((kl*m_eff2)/(k2*m_effl)))*exp(-(i*(xl*(kl +k2)))))); 
m21 =((1/2)*((1 -((kl*m_eff2)/(k2*m_effl)))*exp(i*(xl*(kl +k2))))); 
m22 = ((1/2)*((1 + ((kl*m_eff2)/(k2*m_effl)))*exp(-(i*(xl*(kl - k2)))))); 

m_ll = ((1/2)*((1 + ((k2*m_eff3)/(k3*m_eff2)))*exp(i*(x2*(k2 - k3))))); 
m_12 = ((1/2)*((1 - ((k2*m_eff3)/(k3*m_eff2)))*exp(-(i*(x2*(k2 + k3)))))); 
m_21 = ((1/2)*((1 - ((k2*m_efB)/(k3*m_eff2)))*exp(i*(x2*(k2 +k3))))); 
m_22 = ((1/2)*((1 + ((k2*m_efB)/(k3*m_eff2)))*exp(-(i*(x2*(k2 - k3)))))); 

m_ll_ = ((1/2)*((1 + ((k3*m_eff4)/(k4*m_efB)))*exp(i*(x3*(k3 - k4))))); 
m_12_ = ((1/2)*((1 - ((k3*m_eff4)/(k4*m_efB)))*exp(-(i*(x3*(k3 + k4)))))); 
m_21_ = ((1/2)*((1 - ((k3*m_eff4)/(k4*m_efB)))*exp(i*(x3*(k3 + k4))))); 
m_22_ = ((1/2)*((1 + ((k3*m_eff4)/(k4*m_efB)))*exp(-(i*(x3*(k3 - k4)))))); 


Ml =[mll,ml2;m21,m22]; 

M2 = [m_ll,m_12;m_21,m_22]; 

M3 = [m_l l_,m_12_;m_21_,m_22_]; 


%Now we find an expression for M22 


B = M3*M2*M1; 
J22 = B(2,2); 


%Program used to calculate and plot the wavefunctions 


Psiaqw 

%Defme the function y- 


z = y; 

y = real(J22); 

%Define condition for energy eigenvalues- 

ifz.*y <0 
disp(E) 

Eout(a) = E; %Stores energy levels into an array 
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a = a+1; 

%plottmg the wavefunctions- 

hold on 

plot((X*lel0),((psil/20)+E)) 
plot((Y*lel0),((psi2/20)+E)) 
plot((Z*lel0),((psi3/20)+E)) 
plot((T* 1 e 10),((psi4/20)+E)) 

end 

end 

%This program finds the energy levels in the Valence Band 
valenceInGaAsP 

%Drawing the finite well- 

well2 %This program draws the Conduction band well 

%-End of aqw InGaAsP- 
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%This program is titled Psi aqw 
%Created: January 2002 
%Author: Kevin Lantz 


%This is a program that is used to find the wavefunctions in each region of the quantum well. 
%The bulk of this program is trying to find the normalization constant hi, which determines the 
%normalized wavefunctions and the coefficients A1,B1,A2,B2,A3,B3,A4,B4... 

deltax = le-10; 

%A11 of the coefficients can be written in terms of B1 

B1 = 1; %dummy variable...actual B1 is Bl*bl == hi 
A2 = ml2*Bl; 

B2 = m22*Bl; 

A3 = (m_ll*ml2 + m_12*m22)*Bl; 

B3 = (m_21*ml2 + m_22*m22)*Bl; 

A4 = (m_ll_*(m_ll*ml2 + m_12*m22) + m_12_*(m_21*ml2 + m_22*m22))*Bl; 

%Defming the regions of the quantum well- 

X = -(50*deltax):deltax:0; 

Y = 0:deltax:L2; 

Z = L2:deltax:L3; 

T = L3:deltax:L3+(50*deltax); 

%Defining the unnormalized wavefunctions- 

psil = Bl*exp(-(i*kl*X)); 
psi2 = A2*exp(i*k2*Y) + B2*exp(-(i*k2*Y)); 
psi3 = A3*exp(i*k3*Z) + B3*exp(-(i*k3*Z)); 
psi4 = A4*exp(i*k4*T); 

%Integrating- 

psil = conj(psil); 
psi_2 = conj(psi2); 
psi_3 = conj(psi3); 
psi_4 = conj(psi4); 

% the integral over all space of |psi|^2 == psi*conj(psi) has to be 1 

integrandl = (psil.*psi_l)*deltax; 
integrand2 = (psi2.*psi_2)*deltax; 
integrand3 = (psi3.*psi_3)*deltax; 
integrand4 = (psi4.*psi_4)*deltax; 

integral 1 = sum(integrandl); 
integral2 = sum(integrand2); 
integraB = sum(integrand3); 
integral4 = sum(integrand4); 

fmali = integral 1 + integraB + integraB + integraM; 

hi = l/(sqrt(fmali)); %normahzation constant 

% normalized wavefunctions are just the unnormalized wavefunctions multiplied by the 
normalization constant 
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Psil = bl*psil; 

Psi2 = bl*psi2; 

Psi3 = bl*psi3; 

Psi4 = bl*psi4; 

% This part of the program integrates the normalized wavefunction in order to check 
%that it equals 1. 

%This is used to make sure that the program is working... 

Psil = conj(Psil); 

Psi_2 = conj(Psi2); 

Psi_3 = conj(Psi3); 

Psi_4 = conj(Psi4); 

Integrandl = (Psil.*Psi_l)*deltax; 

Integrand2 = (Psi2.*Psi_2)*deltax; 

Integrand3 = (Psi3.*Psi_3)*deltax; 

Integrand4 = (Psi4.*Psi_4)*deltax; 

Integral 1 = sum(Integrandl); 

Integral2 = sum(Integrand2); 

IntegraB = sum(Integrand3); 

IntegraM = sum(Integrand4); 

fmall = Integral 1 + IntegraB + IntegraB + IntegraM; 

%-End of Psi aqw- 
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%This program is titled E InGaAsP 
%Created: February 2002 
%Author: Kevin Lantz 

%This program has the same form as aqw InGaAsP except it uses Airy functions 
clc 

F = input('Electric Field in [V/m]: '); 

Xl_ = input('Percentage of A1 (0-1) in region 1: '); 

X2_ = input('Percentage of A1 (0-1) in region 2: '); 

X3_ = input('Percentage of A1 (0-1) in region 3: '); 

X4_ = X1_; 

width 1 = Inf; %Extends from -Infinity to 0; 

width2 = input('Width of the region 0->a in angstroms: '); %Well width 

width3 = input('Width of the region a->b in angstroms: '); 

width4 = Inf; %Extends from L to Infinity; 


% - 

y = 0; %InitiaI value of y to correct the logic of the later condition 

a= 1; 

% - % 


El = widthl*le-10; %Change input of angstroms->meters 
L2 = width2*le-10; 

L3 = (width2-t-width3)*le-10; 

L4 = width4*le-10; 

hbar= 1.055e-34; %[J*s] 
m0 = 9.11e-31;%[kg] 
e= 1.602e-19;%[C] 

%Conditions used to make the program function correctly under positive or negative bias 
ifF<0; 
zl =L3; 
z2 = L2; 
z3 = 0; 

XI =X4_; 

X2 = X3_; 

X3 = X2_; 

X4 = X1_; 
else 
zl =0; 
z2 = L2; 
z3 = L3; 

XI =X1_; 

X2 = X2_; 

X3 = X3_; 

X4 = X4_; 
end 

% - % 

ifF<0 

egl = 1.4236*(1-X1) 0.81612*X1 - 0.13*(X1*(1-X1)); 

eg3 = 1.519*(1-X2) 0.417*X2 - 0.477*(X2*(1-X2)); 

eg2 = 1.4236*(1-X3) 0.81612*X3 - 0.13*(X3*(1-X3)); 
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eg4 = 1.4236*(1-X1) + 0.81612*X1 - 0.13*(X1*(1-X1)); 
m_effl = l/((Xl/(0.0365*m0)) + ((l-Xl)/(0.0795*m0))); 
m_efG = l/((X2/(0.026*m0)) + ((l-X2)/(0.067*m0))); 
mem = l/((X3/(0.0365*m0)) + ((l-X3)/(0.0795*m0))); 
m_eff4 = l/((Xl/(0.0365*m0)) + ((l-Xl)/(0.0795*m0))); 


else 

egl = 1.4236*(1-X1) + 0.81612*X1 - 0.13*(X1*(1-X1)); 
eg2 = 1.519*(1-X2) + 0.417*X2 - 0.477*(X2*(1-X2)); 
eg3 = 1.4236*(1-X3) + 0.81612*X3 - 0.13*(X3*(1-X3)); 
eg4 = 1.4236*(1-X1) + 0.81612*X1 - 0.13*(X1*(1-X1)); 
m_effl = l/((Xl/(0.0365*m0)) + ((l-Xl)/(0.0795*m0))); 
m_eff2 = l/((X2/(0.026*m0)) + ((l-X2)/(0.067*m0))); 
m_eff3 = l/((X3/(0.0365*m0)) + ((l-X3)/(0.0795*m0))); 
m_eff4 = l/((Xl/(0.0365*m0)) + ((l-Xl)/(0.0795*m0))); 
end 


VI =0.6*(egl -0.8161193); 

V2 = 0.6*(eg2 - 0.8161193); %Clioose the origin as the bottom left comer of region 2. 
V3 = 0.6*(eg3 -0.8161193); 

V4 = V1; 

% - % 

Cl = (((2*m_effl*e)/((F*hbar)^2))^(l/3)); 

C2 = (((2*m_eff2*ey((F*hbar)^2))^(l/3)); 

C3 = (((2*m_efG*e)/((F*hbar)^2))^(l/3)); 

C4 = (((2*m_eff4*e)/((F*hbar)^2))^(l/3)); 

sigmal = ((m_effl*C2)/(m_eff2*Cl)); 
sigma2 = ((m_eff2*C3)/(m_efG*C2)); 
sigma3 = ((m_efG*C4)/(m_eff4*C3)); 


% - % 

ifF<0; 

DD = 0; 

V2_ = V1; 
else 

DD = -F*L2; 

V2_ = V1 -F*L3; 
end 

% - % 

forE=DD:0.001:V2_; 

etal = (E - VI); 
eta2 = (E - V2); 
eta3 = (E - V3); 
eta4 = (E - V4); 


Alphal = -Cl*(F*zl -I- etal); 
Alpha2 = -C2*(F*zl -I- eta2); 
Alpha2_ = -C2*(F*z2 -I- eta2); 
Alpha3 = -C3*(F*z2 -I- eta3); 
Alpha3_ = -C3*(F*z3 -I- eta3); 
Alpha4 = -C4*(F*z3 -I- eta4); 


%■ 


% 
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Ewell2 

o/o- 


mil = (pi)*(airy(Alpha2)*airy(3,Alphal) - sigmal*(airy(2,Alphal)*airy(l,Alpha2))); 
ml2 = (pi)*(airy(2,Alpha2)*airy(3,Alphal) - sigmal*(airy(2,Alphal)*airy(3,Alpha2))); 
m21 = (pi)*(sigmal*(airy(Alphal)*airy(l,Alpha2)) - airy(l,Alphal)*airy(Alpha2)); 
m22 = (pi)*(sigmal*(airy(Alphal)*airy(3,Alpha2)) - airy(l,Alphal)*airy(2,Alpha2)); 

m il = (pi)*(airy(Alpha3)*airy(3,Alpha2_) - sigma2*(airy(2,Alpha2_)*airy(l,Alpha3))); 
m_12 = (pi)*(airy(2,Alpha3)*airy(3,Alpha2_) - sigma2*(airy(2,Alpha2_)*airy(3,Alpha3))); 
m_21 = (pi)*(sigma2*(airy(Alpha2_)*airy(l,Alpha3)) - airy(l,Alpha2_)*airy(Alpha3)); 
m_22 = (pi)*(sigma2*(airy(Alpha2_)*airy(3,Alpha3)) - airy(l,Alpha2_)*airy(2,Alpha3)); 

m_ll_ = (pi)*(airy(Alpha4)*airy(3,Alpha3_) - sigma3*(airy(2,Alpha3_)*airy(l,Alpha4))); 
m_12_ = (pi)*(airy(2,Alpha4)*airy(3,Alpha3_) - sigma3*(airy(2,Alpha3_)*airy(3,Alpha4))); 
m_21_ = (pi)*(sigma3*(airy(Alpha3_)*airy(l,Alpha4)) - airy(l,Alpha3_)*airy(Alpha4)); 
m_22_ = (pi)*(sigma3*(airy(Alpha3_)*airy(3,Alpha4)) - airy(l,Alpha3_)*airy(2,Alpha4)); 

Ml =[mll,ml2;m21,m22]; 

M2 = [m_l l,m_12;m_21,m_22]; 

M3 = [m_ll_,m_12_;m_21_,m_22J; 

M2M3 = M2*M3; 

M = M1*M2*M3; 

M22 = real(M(2,2)); 

Epsiaqw 

z = y; 
y = M22; 

ifz.*y <0 
disp(E) 

Eout(a) = E; 
a = a + 1; 
hold on 

plot((X* 1 e 10),((Psi 1 ./(4e5))+E)) 
plot((Y*lelO),((Psi2./(4e5))+E)) 
plot((Z * 1 e 10), ((Psi3 ./(4e5))+E)) 
plot((T* 1 e 10),((Psi4./(4e5))+E)) 


end 

end 

E valence InGaAsP 


■End of E InGaAsP- 
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%This program is titled absorption. 

%Created: February 2002 
%Author: Kevin Lantz 

%This program uses the wavefunctions to find the bound-to-bound and bound-to-continuum 
%absorption spectra. 

clc 

aqwIuGaAsP 


%-Psi at Eout(l) 


kEl_l = sqrt((2*m_effl*(Eout(l)-Vl)*q)/(hbar^2)); 
kEl_2 = sqrt((2*m_eff2*(Eout(l)-V2)*q)/(hbar^2)); 
kEl_3 = sqrt((2*m_eff3*(Eout(l)-V3)*q)/(hbar^2)); 
kEl_4 = sqrt((2*m_eff4*(Eout(l)-V4)*q)/(hbar^2)); 

mEl_ll = ((1/2)*((1 -I- ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(i*(xl*(kEl_l -kEl_2))))); 
mEl_12 = ((1/2)*((1 - ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(-(i*(xl*(kEl_l -I- kEl_2)))))); 
mEl_21 = ((1/2)*((1 - ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(i*(xl*(kEl_l -I- kEl_2))))); 
mEl_22 = ((1/2)*((1 -I- ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(-(i*(xl*(kEl_l - kEl_2)))))); 

m_El_ll = ((1/2)*((1 -I- ((kEl_2*m_efB)/(kEl_3*m_eff2)))*exp(i*(x2*(kEl_2 - kEl_3))))); 
m_El_12 = ((1/2)*((1 - ((kEl_2*m_efB)/(kEl_3*m_eff2)))*exp(-(i*(x2*(kEl_2 -I- kEl_3)))))); 
m_El_21 = ((1/2)*((1 - ((kEl_2*m_efB)/(kEl_3*m_eff2)))*exp(i*(x2*(kEl_2 -I- kEl_3))))); 
m_El_22 = ((1/2)*((1 -I- ((kEl_2*m_eff3)/(kEl_3*m_eff2)))*exp(-(i*(x2*(kEl_2 - kEl_3)))))); 

m_El_l 1_ = ((1/2)*((1 -I- ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(i*(x3*(kEl_3 - kEl_4))))); 
m_El_12_ = ((1/2)*((1 - ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(-(i*(x3*(kEl_3 -I- kEl_4)))))); 
m_El_21_ = ((1/2)*((1 - ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(i*(x3*(kEl_3 -I- kEl_4))))); 
m_El_22_ = ((1/2)*((1 -I- ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(-(i*(x3*(kEl_3 - kEl_4)))))); 

ME1_1 = [mEl_ll,mEl_12;mEl_21,mEl_22]; 

ME1_2 = [m_El_ll,m_El_12;m_El_21,m_El_22]; 

ME1_3 = [m_El_ll_,m_El_12_;m_El_21_,m_El_22J; 

M2M1 =ME1_2*ME1_1; 

M3M2M1 =ME1_3*ME1_2*ME1_1; 

AE1_1 = 0; 

BE1_1 = 1; 

AE1_2 = ME1_1(1,1)*AE1_1 ME1_1(1,2)*BE1_1; 

BE1_2 = ME1_1(2,1)*AE1_1 ME1_1(2,2)*BE1_1; 

AE1_3 = M2M1(1,1)*AE1_1 M2M1(1,2)*BE1_1; 

BE1_3 =M2M1(2,1)*AE1_1 -HVI2M1(2,2)*BE1_1; 

AE1_4 = M3M2M1(1,1)*AE1_1 M3M2M1(1,2)*BE1_1; 

BE1_4 = 0; 


deltax = le-10; 


X = -L3:deltax:0; 

Y = 0:deltax:L2; 

Z = L2:deltax:L3; 

T = L3:deltax:(2*L3); 

psiEl l = AEl_l*exp(i*kEl_l*X) -I- BEl_l*exp(-(i*kEl_l*X)); 
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psiEl_2 = AEl_2*exp(i*kEl_2*Y) + BEl_2*exp(-(i*kEl_2*Y)); 
psiEl_3 = AEl_3*exp(i*kEl_3*Z) + BEl_3*exp(-(i*kEl_3*Z)); 
psiEl_4 = AEl_4*exp(i*kEl_4*T) + BEl_4*exp(-(i*kEl_4*T)); 

psiEll = conj(psiEl_l); 
psi_El_2 = conj(psiEl_2); 
psi_El_3 = conj(psiEl_3); 
psi_El_4 = conj(psiEl_4); 

deltalEll = (psi_El_l.*psiEl_l)*deltax; 
deltaIEl_2 = (psi_El_2.*psiEl_2)*deltax; 
deltaIEl_3 = (psi_El_3.*psiEl_3)*deltax; 
deltaIEl_4 = (psi_El_4.*psiEl_4)*deltax; 

IE 11 = sum(deltaIEl_l); 

IE1_2 = sum(deltaIEl_2); 

IE1_3 = sum(deltaIEl_3); 

IE1_4 = sum(deltaIEl_4); 

finalEl = IE1_1+IE1_2+IE1_3+IE1_4; 

bEl = l/(sqrt(finalEl)); 

PsiEl_l =bEl*psiEl_l; 

PsiEl_2 = bEl*psiEl_2; 

PsiEl_3=bEl*psiEl_3; 

PsiEl_4 = bEl*psiEl_4; 

%-derivative of Psi at Eout(l)- 

psiEEl l = AEl_l*exp(i*kEl_l*X) - BEl_l*exp(-(i*kEl_l*X)); 
psiEEl_2 = AEl_2*exp(i*kEl_2*Y) - BEl_2*expHi*kEl_2*Y)); 
psiEEl_3 = AEl_3*exp(i*kEl_3*Z) - BEl_3*exp(-(i*kEl_3*Z)); 
psiEEl_4 = AEl_4*exp(i*kEl_4*T) - BEl_4*exp(-(i*kEl_4*T)); 

PsiEEll =bEl*psiEEl_l; 

PsiEEl_2 = bEl*psiEEl_2; 

PsiEEl_3 =bEl*psiEEl_3; 

PsiEEl_4 = bEl*psiEEl_4; 

%-Psi at Eout(2)- 

kE2_l = sqrt((2*m_effl*(Eout(2)-Vl)*q)/(hbari'2)); 
kE2_2 = sqrt((2*m_eff2*(Eout(2)-V2)*q)/(hbar^2)); 
kE2_3 = sqrt((2*m_efO*(Eout(2)-V3)*q)/(hbari'2)); 
kE2_4 = sqrt((2*m_eff4*(Eout(2)-V4)*q)/(hbar^2)); 

mE2_ll = ((1/2)*((1 + ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(i*(xl*(kE2_l - kE2_2))))); 
mE2_12 = ((1/2)*((1 - ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(-(i*(xl*(kE2_l + kE2_2)))))); 
mE2_21 = ((1/2)*((1 - ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(i*(xl*(kE2_l +kE2_2))))); 
mE2_22 = ((1/2)*((1 + ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(-(i*(xl*(kE2_l - kE2_2)))))); 

m_E2_ll = ((1/2)*((1 + ((kE2_2*m_eff3)/(kE2_3*m_eff2)))*exp(i*(x2*(kE2_2 - kE2_3))))); 
m_E2_12 = ((1/2)*((1 - ((kE2_2*m_efG)/(kE2_3*m_eff2)))*exp(-(i*(x2*(kE2_2 + kE2_3)))))); 
m_E2_21 = ((1/2)*((1 - ((kE2_2*m_efG)/(kE2_3*m_eff2)))*exp(i*(x2*(kE2_2 +kE2_3))))); 
m_E2_22 = ((1/2)*((1 + ((kE2_2*m_eff3)/(kE2_3*m_eff2)))*exp(-(i*(x2*(kE2_2 - kE2_3)))))); 
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m_E2_l 1_ = ((1/2)*((1 + ((kE2_3*m_eff4)/(kE2_4*m_efD)))*exp(i*(x3*(kE2_3 - kE2_4))))); 
m_E2_12_ = ((1/2)*((1 - ((kE2_3*m_eff4)/(kE2_4*m_efG)))*exp(-(i*(x3*(kE2_3 + kE2_4)))))); 
m_E2_21_ = ((1/2)*((1 - ((kE2_3*m_eff4)/(kE2_4*m_efG)))*exp(i*(x3*(kE2_3 +kE2_4))))); 
m_E2_22_ = ((1/2)*((1 + ((kE2_3*m_eff4)/(kE2_4*m_efG)))*exp(-(i*(x3*(kE2_3 - kE2_4)))))); 


ME2_1 = [mE2_ll,mE2_12;mE2_21,mE2_22]; 

ME2_2 = [m_E2_ll,m_E2_12;m_E2_21,m_E2_22]; 

ME2_3 = [m_E2_ll_,m_E2_12_;m_E2_21_,m_E2_22J; 

M2_M1 =ME2_2*ME2_1; 

M3M2_M1 =ME2_3*ME2_2*ME2_1; 

AE2_1 = 0; 

BE2_1 = 1; 

AE2_2 = ME2_1(1,1)*AE2_1 + ME2_1(1,2)*BE2_1; 

BE2_2 = ME2_1(2,1)*AE2_1 + ME2_1(2,2)*BE2_1; 

AE2_3 = M2_M1(1,1)*AE2_1 + M2_M1(1,2)*BE2_1; 

BE2_3 = M2_M1(2,1)*AE2_1 + M2_M1(2,2)*BE2_1; 

AE2_4 = M3M2_M1(1,1)*AE2_1 + M3M2_M1(1,2)*BE2_1; 
BE2_4 = 0; 

psiE2_l = AE2_l*exp(i*kE2_l*X) + BE2_l*exp(-(i*kE2_l*X)); 
psiE2_2 = AE2_2*exp(i*kE2_2*Y) + BE2_2*expHi*kE2_2*Y)); 
psiE2_3 = AE2_3*exp(i*kE2_3*Z) + BE2_3*exp(-(i*kE2_3*Z)); 
psiE2_4 = AE2_4*exp(i*kE2_4*T) + BE2_4*expHi*kE2_4*T)); 


psi_E2_l = conj(psiE2_l); 
psi_E2_2 = conj(psiE2_2); 
psi_E2_3 = conj(psiE2_3); 
psi_E2_4 = conj(psiE2_4); 


deltaIE2_l = (psi_E2_l.*psiE2_l)*deltax; 
deltaIE2_2 = (psi_E2_2.*psiE2_2)*deltax; 
deltaIE2_3 = (psi_E2_3.*psiE2_3)*deltax; 
deltaIE2_4 = (psi_E2_4.*psiE2_4)*deltax; 


IE2_1 = sum(deltaIE2_l); 

IE2_2 = sum(deltaIE2_2); 

IE2_3 = sum(deltaIE2_3); 

IE2_4 = sum(deltaIE2_4); 

finalE2 = IE2_1+IE2_2+IE2_3+IE2_4; 

bE2 = l/(sqrt(finalE2)); 

PsiE2_l =bE2*psiE2_l; 

PsiE2_2 = bE2*psiE2_2; 
PsiE2_3=bE2*psiE2_3; 

PsiE2_4 = bE2*psiE2_4; 

PsiE2_l_ = conj(PsiE2_l); 

PsiE2_2_ = conj(PsiE2_2); 

PsiE2_3_ = conj(PsiE2_3); 

PsiE2_4_ = conj(PsiE2_4); 
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%■ 


■<Psi(Eout(2))*|k|Psi(Eout(l))>- 


deltaE12_l = (PsiE2_l_.*kEl_l.*PsiEEl_l)*deltax; 
deltaE12_2 = (PsiE2_2_.*kEl_2.*PsiEEl_2)*deltax; 
deltaE12_3 = (PsiE2_3_.*kEl_3.*PsiEEl_3)*deltax; 
deltaE12_4 = (PsiE2_4_.*kEl_4.*PsiEEl_4)*deltax; 

E12_l = sum(deltaE12_l); 

E12_2 = sum(deltaE12_2); 

E12_3 = sum(deltaE12_3); 

E12_4 = sum(deltaE12_4); 

averageE12 = E12_l+E12_2+E12_3+E12_4; 

AverageE12 = averageE12*(conj(averageE12)); 

%-Line Shape- 

N= le24; 
eO = 8.85e-12; 
c = 3e8; 
n = 3.5; 

Gamma = ((10e-3)*q)/2; 
ap= 1; 

for EE = (Eout(l)-0.1):le-4:(Vl-0.001); 

g = (Gamma)/((((Eout(2)-Eout(l)-EE)*q)^2)+((Gamma)^2)); 
gmax = (Gamma)/((((0)*q)^2)+((Gamma)^2)); 

C = (N*(q^2)*(hbar^3))/((m_effl^2)*e0*n*c); 
alphaE12(ap) = (C.*g.*AverageE12)/(EE*q); 

ap = ap + 1; 
end 


% 

%-Psi above the well- 

r=l; 

for Energy = (Vl+0.0001):0.0001:(Vl+.025) 


kAl = sqrt((2*m_effl*(Energy-Vl)*q)/(hbar''2)); 
kA2 = sqrt((2*m_eff2*(Energy-V2)*q)/(hbar''2)); 
kA3 = sqrt((2*m_eff3*(Energy-V3)*q)/(hbar''2)); 
kA4 = sqrt((2*m_eff4*(Energy-V4)*q)/(hbar''2)); 


mA_ll = ((1/2)*((1 + ((kAl*m_eff2)/(kA2*m_effl)))*exp(i*(xl*(kAl - kA2))))); 
mA_12 = ((1/2)*((1 - ((kAl*m_eff2)/(kA2*m_effl)))*exp(-(i*(xl*(kAl +kA2)))))); 
mA_21 = ((1/2)*((1 - ((kAl*m_eff2)/(kA2*m_effl)))*exp(i*(xl*(kAl +kA2))))); 
mA_22 = ((1/2)*((1 + ((kAl*m_eff2)/(kA2*m_effl)))*exp(-(i*(xl*(kAl - kA2)))))); 

m_A_ll = ((1/2)*((1 + ((kA2*m_efG)/(kA3*m_eff2)))*exp(i*(x2*(kA2 - kA3))))); 
m_A_12 = ((1/2)*((1 - ((kA2*m_efG)/(kA3*m_eff2)))*exp(-(i*(x2*(kA2 +kA3)))))); 
m_A_21 = ((1/2)*((1 - ((kA2*m_efG)/(kA3*m_eff2)))*exp(i*(x2*(kA2 +kA3))))); 
m_A_22 = ((1/2)*((1 + ((kA2*m_efG)/(kA3*m_eff2)))*exp(-(i*(x2*(kA2 - kA3)))))); 

m_A_ll_ = ((1/2)*((1 + ((kA3*m_eff4)/(kA4*m_efB)))*exp(i*(x3*(kA3 - kA4))))); 

63 









m_A_12_ = ((1/2)*((1 - ((kA3*m_eff4)/(kA4*m_efB)))*exp(-(i*(x3*(kA3 +kA4)))))); 
m_A_21_ = ((1/2)*((1 - ((kA3*m_eff4)/(kA4*m_efB)))*exp(i*(x3*(kA3 +kA4))))); 
m_A_22_ = ((1/2)*((1 + ((kA3*m_eff4)/(kA4*m_efB)))*exp(-(i*(x3*(kA3 - kA4)))))); 

MA_1 = [mA_ll,mA_12;mA_21,mA_22]; 

MA_2 = [m_A_ll,m_A_12;m_A_21,m_A_22]; 

MA_3 = [m_A_ll_,m_A_12_;m_A_21_,m_A_22J; 


M2M1_ = MA_2*MA_1; 

M3M2M1_ = MA_3*MA_2*MA_1; 

A1 = 1; 

B1 = ((-M3M2M1_(2,1))/(M3M2M1_(2,2)))*A1; 
A2 = MA_1(1,1)*A1 + MA_1(1,2)*B1; 

B2 = MA_1(2,1)*A1 + MA_1(2,2)*B1; 

A3 = M2M1_(1,1)*A1 + M2M1_(1,2)*B1; 

B3 = M2M1_(2,1)*A1 + M2M1_(2,2)*B1; 

A4 = M3M2M1_(1,1)*A1 + M3M2M1_(1,2)*B1; 
B4 = 0; 


PsiAl = Al*exp(i*kAl*X) + Bl*exp(-(i*kAl*X)); 
PsiA2 = A2*exp(i*kA2*Y) + B2*exp(-(i*kA2*Y)); 
PsiA3 = A3*exp(i*kA3*Z) + B3*exp(-(i*kA3*Z)); 
PsiA4 = A4*exp(i*kA4*T) + B4*exp(-(i*kA4*T)); 

PsiAl_ = conj(PsiAl); 

PsiA2_ = conj(PsiA2); 

PsiA3_ = conj(PsiA3); 

PsiA4_ = conj(PsiA4); 

%-<Psi(above well)* |k|Psi(Eout( 1 ))>- 


deltaEAl = (PsiAl_.*kEl_l.*PsiEEl_l)*deltax; 
deltaEA_2 = (PsiA2_.*kEl_2.*PsiEEl_2)*deltax; 
deltaEA_3 = (PsiA3_.*kEl_3.*PsiEEl_3)*deltax; 
deltaEA_4 = (PsiA4_.*kEl_4.*PsiEEl_4)*deltax; 

EA_1 = sum(deltaEA_l); 

EA_2 = sum(deltaEA_2); 

EA_3 = sum(deltaEA_3); 

EA_4 = sum(deltaEA_4); 

averageEA = EA_1+EA_2+EA_3+EA_4; 

AverageEA(r) = averageEA* (conj (averageEA)); 

%-Finding the absorption coefficient- 

C = (N*(q^2)*(hbar^2))/(2*(m_effB2)*eO*n*c); 

CC = (((((Energy-Vl))*q)^2) + (Gamma^2)); 

rho = ((l/CC)^(l/4)); 

C_ = (((Energy-Eout(l))*q)^2); 
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A_ = (sqrt((2*m_effl)/C_))*rho; 

alpha(r) = C*A_*AverageEA(r); 

r = r+ 1; 
end 

%-End of absorption- 
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